A twist localizes three-dimensional patterns 
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A mechanism for the localization of spatially periodic, self-organized patterns in anisotropic media 
which requires systems extended in all three spatial dimensions is presented: When the anisotropy 
axis is twisted the pattern becomes localized in planes parallel to the anisotropy axis. An analytic 
description of the effect is developed and used to interpret recent experiments in the high-frequency 
regime of electroconvection by Bohatsch and Stannarius [Phys. Rev. E 60, 5591 (1999)]. The 
localization width is found to be of the order of magnitude of the geometrical average of pattern 
wavelength and the inverse twist. 
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I. INTRODUCTION 



Three-dimensional pattern formation in dissipative 
systems has been theoretically investigated to some ex- 
tent for the case of the complex Ginzburg-Landau equa- 
tion — but there seems to be no experimental 
model for which the theory is quantitatively valid. Qual- 
itatively, many results should also apply to a similar 
case, pattern formation in 3D excitable media, for which 
controlled experimental realizations exist (e.g. @-||]), al- 
though the containers are typically rather small. In both 
systems it seems reasonable to concentrate on the struc- 
ture and dynamics of the filaments formed by the cores 
of scroll waves, which largely control the dynamics in 
the remaining space. In a recent theoretical work, for 
example, the filament dynamics was investigated for an 
anisotropic excitable medium with a twisted anisotropy 
axis pO[ , a situation found in the human heart muscle. 
It was predicted that the filament typically drifts such as 
to align with the anisotropy axis. 

Both the complex Ginzburg-Landau equation and ex- 
citable media exhibit traveling waves. For systems form- 
ing stationarily bifurcating spatially periodic structures 
the dominance of the filaments formed by topological 
defects is not so strong, and different effects can be ex- 
pected. A prototype for such a system might be the high- 
frequency regime of electroconvection (EC) in nematic 
liquid crystals. Convection patterns with a wavelength 
A which is small compared to the spatial extensions of 
the sample are easily obtained and sustained over long 
times by applying a sufficiently high ac voltage (with a 
modulation frequency higher than the "cutoff frequency" 
« the inverse charge relaxation time) at plane, parallel 
(and transparent) electrodes enclosing the liquid crystal. 
However, there is still some discussion on whether the 
patterns are really three dimensional, i.e. arising from 
a bulk instability and forming a laminar structure, or 



rather parallel rolls located at the electrodes. 

Closely related to this geometrical question is the ques- 
tion about the mechanism of convection. The most pop- 
ular explanation, here called the standard model (SM), 
is by the dielectric regime of EC via the Carr-Helfrich 
mechanism (for details, see the references below), which 
is based the Leslie-Ericksen formulation of nematody- 
namics. The SM predicts bulk convection. It has, on 
the linear level, been worked out by Dubois- Violette et 
al. [ pT| . Later it was modified for square- wave dri ving 
voltages and to include flexoelectric effects p^-[l5|] 
and the effects of finite cell thickness [16[. Recently it 
was extended to the nonlinear level ||l7[~^But attempts 
to bring the predictions of the SM (threshold voltage and 
critical wave number) into quantitative agreement with 
experiments were only partly successful. 

The main challenge to the SM is the "isotropic" mech- 
anism for EC ||l8|,|9|, for which the anisotropy of the 
nematic is only of secondary relevance. There, EC is 
driven by inverted charge gradients in thin Debye layers 
near the electrodes, which are not included in the SM. 
The "isotropic" mechanism naturally leads to boundary 
convection. Since several material parameters relevant to 
the "isotropic" mechanism are unknown, only semiquan- 
titative predictions of the effect can be made. 

In principle, both mechanisms are active, but it is rea- 
sonable to assume that one of them dominates. In a re- 
cent effort to settle the controversy on bulk vs. boundary 
instability, Bohatsch and Stannarius |^| developed an 
experiment with a modified design, which was intended 
to map the structure across the cell onto the 2D shadow- 
graph image (or the diffraction pattern) of light passing 
through the cell. 

The technique of Bohatsch and Stannarius was to used 
samples with a twisted geometry: just as in TN LCD dis- 
plays, the anchoring of the nematic director fi (n^ — 1) 
at the upper and lower electrodes (separation d) is such 
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that n is parallel to the layer but the alignment axis of the 
director at the upper boundary is perpendicular to that 
at the lower boundary. In the equilibrium configuration 
the director is planarly oriented and interpolates between 
the boundaries with a uniform twist n ■ rot fi. Configura- 
tions with positive or negative twist are thinkable. They 
are related to each other by a reflection. Without loss 
of generality a negative twist {—it/ 2d) shall be assumed. 
Although inversion symmetry is broken in the twisted 
geometry, there is still the symmetry with respect to a 
180° rotation of the cell around the axis of the direc- 
tor at mid plane. The critical mode at the onset of EC, 
which breaks the approximate translation symmetry, can 
be symmetric with respect to this rotation (perhaps fol- 
lowed by a translation by half a wavelength), leading to 
normal rolls, or break also this discrete symmetry, then 
leading to oblique rolls. 

In experiments with untwisted cells, the rolls (or lam- 
inae) observed at the threshold of structure formation in 
the high-frequency regime are always normal^ ^. Thus 
Bohatsch and Stannarius conclude by an intuitive argu- 
ment that, in the twisted cells, the observation of nor- 
mal rolls at threshold demonstrates a bulk instability, 
while the observation of oblique rolls is an indicator for a 
boundary instability, in particular when the rolls are ap- 
proximately normal to the alignment of n at the bound- 
aries. They observe a transition from normal to oblique 
rolls as the ac frequency increases. 

Assuming that the SM is basically correct. Sections || 
and ^ of this paper develop the linear theory for the 
convection pattern in twisted cells, predicting critical 
modes localized in planes far from or near to the elec- 
trodes, respectively, depending on the boundary condi- 
tions. Section [V sketches the nonlinear effects relevant 
at the threshold of convection and Section ^ concludes, 
in particular with the interpretation of the experiments 
in the light of the theory. 



II. LOCALIZED BULK MODES IN THE 
STANDARD MODEL 



Below the linear theory of high-frequency EC in 
twisted cells is developed within the framework of the 
SM. The association of the obliqueness angle with the po- 
sition of the convective mode turns out to be correct. But 
it is also found that the experimental results for twisted 
cells cannot be transfered to the untwisted cell as it was 
suggested by Bohatsch and Stannarius, because the lo- 
calization of the linear mode, near the boundaries or in 



the bulk, is enforced by the twist itself. 

The calculations are based on the general weakly non- 
linear theory of dielectric EC in large containers, 
which has to take the interaction of the convection pat- 
tern with several homogeneous soft modes (hydrody- 
namic modes) explicitly into account. One of these 
modes is the director orientation. As follows from the 
general calculation, a twist does, at lowest order, not ex- 
cite the other soft modes of the liquid crystal, such as 
electric potentials or fluid motion, even when the electric 
ac voltage is applied. The direct influence of the twist 
on the convection pattern reduces to a purely geomet- 
ric effect (no additional coupling constants not relevant 
for the untwisted cell are involved). The local amplitude 
and phase of the convection pattern is expressed by a 
complex field A — A{x, y, z, t) such that the modulation 
of actual physical quantities u is given by u{x,y,z,t) = 
^0^(2;, y, z, t) exp(i5c2^)+c.c.-|-(higher order corrections), 
where qc is the wave number of the convection pattern 
at threshold in the untwisted cell and uq is a complex 
constant. The linear dynamics of A in the presence of 
a pure, not too strong twist of the nematic director is 
described by 

rdtA = [e + ildl + ^{dy - iq.ipf + e.dl] A. (1) 

Here the z axis is perpendicular to the electrodes. To 
be specific assume that z points "up" and —z points 
"down" . For now, the orientation of x and y axis rel- 
ative to the director anchoring shall be left unspecified. 
This is possible because of the rotational symmetry un- 
derlying Eq. (|l|) The angle between n and the x-z 
plane is denoted by ^p. Quantitative validity of Eq. (|l|) 
requires that be small and that modulations of A and 

in space and time are slow. In particular Eq. (|l|) is 
never valid across the whole cell, since ip varies over an 
angle 0(1), a problem that will be discussed below. The 
(positive) coherence lengths ^j,, ^y, and are typically 
of the order 0{q~^), and r is of the order of the charge 
relaxation time. The strength of the external stress is ex- 
pressed by the reduced control parameter e := E"^ /El — 1 
in terms of the strength of the external electric ac field 
E = V/d and the threshold field for electroconvection Ec 
for the thick [d — > c»), untwisted cell. |e| is assumed to 
be small. (Observe that this definition of e differs from 
the one used in Ref. , which is based on the threshold 
in the twisted cell with finite d, a quantity to be calcu- 
lated below.) Analytic approximations for all quantities 
involved in the description (^ can be found in Ref. . 

The boundary conditions corresponding to Eq. (^ are 

A[x, y, z, y,t) — at the horizontal boundaries. (2) 



'^When chevron patterns are observed, one is not at the 
threshold of structure formation. 

■^In the low- frequency regime of EC they can be normal or 
oblique, however. 
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From the translation invariance in the x-y plane it 
follows that the critical mode is of the form A = 
a{z) ex.p{iQx+iPy) with real Q and P, i.e. with a pattern 
wave vector q = {qc + Q,P)- We will choose the orien- 
tation of the X axis to be parallel to g, in other words 
P = 0. It is easily seen by separation of variables that 
according to Eq. (|l|) the critical mode, i.e. the first mode 
to become unstable as e increases, has (5 = 0. It follows 
that the twist has no influence on the wavelength of the 
critical mode. 

This result is in agreement with the measurements of 
Bohatsch and Stannarius for the cells of intermediate 
thickness. But in their thinnest and in their thickest 
cell they find shorter wavelengths. While the deviation 
in the thin cell can be understood as a breakdown of the 
3D amplitude formalism due to large spatial gradients of 
A and Lp, the reason for the deviation in the thick cell 
is not clear. In agreement with theoretical expectations 
the wavelength decreases with increasing frequency. The 
observed deviation from the approximate A ~ law 
at the highest frequencies indicates that the wavelength 
has reached the Debye length and charge diffusion effects 
need to be taken into account in the theoretical descrip- 
tion, an effect noticed already by Dubois- Violette et al. 
(and typically to be expected at the 10^ fold cutoff 
frequency |p7t). 

Choose the origin of the z axis such that iy9 = at 
2; = 0. Denote by zq the distance of the lower boundary 
from z = (i.e. the lower boundary is at z = — zq). By 
the preceding construction, the amplitude of the critical 
mode is a function of z alone and satisfies 



(3) 



with Lp — ■Kz/2d. The equation has the form of the time- 
independent Schrodinger equation for the quantum har- 
monic oscillator. 

Ignore, for a moment, the horizontal boundaries, and 
impose instead the B.C. 



A is bounded as l^^l — > oo. 
Then the lowest eigenvalue of Eq. (^ is 

(?c S.Z 



e = £r := 



with a critical mode 

A = exp 



2d 



2w^ 



with a width 



2d^ 



1/2 



(4) 



(5) 



(6) 



(7) 



Now assume dqc ^ 1, which is required for the valid- 
ity of Eq. (0). Then = O(dgc)"^ <C 1 and, within 
the range |z| — C(w) covered by the critical mode, 
if = 0{w/d) = 0{dqc)-^/^ < 1. Thus, for |z| = 0(w), 
the solution (||,||) is a physically valid solution of Eq. (|^). 
For |z| ^ w the critical mode is to high accuracy zero 
and thus a solution of the linear problem, also beyond 
the range of validity of Eq. (|l|). Of course any homoge- 
neous B.C., in particular B.C. (^), are satisfied, provided 
the boundaries are sufhciently far away. Equations (||J|) 
therefore give a physically valid threshold and critical 
mode. 

Next, consider the situation that one of two bound- 
aries, e.g. the lower one, is near z = [i.e. zq = 0(iy)], 
the other boundary far away. Then the lower B.C. be- 
comes 



A^Q at z 



(8) 



Equations (|5|,|6|) correspond to the ground state of the 
harmonic oscillator. 



while the upper B.C. remains effectively 

A is bounded as z ^ oo. (9) 

This eigenvalue problem has to be solved numerically 
(after rescaling, it depends only on the parameters e/sc 
and zq/w); here a shooting method is used. The lowest 
eigenvalue £(zo) as a function of zq is shown in Fig. 
(lower curve). As zq decreases to zero and the maxi- 
mum of the critical mode approaches the boundary, the 
Gaussian shape of A{z) is deformed and the threshold 
e(zo) increases monotonically. In particular one finds 
e(0) — Ssc, the eigenvalue of the first excitation and 
lowest antisymmetric mode of the "harmonic oscillator" , 
with the corresponding eigenmode. Excitations near the 
boundaries have higher thresholds than in the bulk and 
are inhibited. 

Thus the theory predicts that in thick enough cells the 
critical mode is localized within a horizontal layer which 
covers only a small fraction = 0{'w/d) — 0{qcd)~^^^ of 
the cell. The orientation of the critical wave vector q 
and the position of the maximum of the critical mode 
within the sample are coupled, such that q is parallel 
to h at the maximum of A(z). The orientation of the 
critical wave vector q lies between the two orientations 
of director anchoring at the boundaries. Small angles of 
order 0{w/d) between q and director anchoring are for- 
bidden, but otherwise the orientation of q with respect 
to the symmetry axis, i.e. the obliqueness angle, is unde- 
termined. One is thus dealing with a highly degenerate 
problem with critical wave vectors lying on a circular arc, 
similar to pattern formation in isotropic systems, where 
the critical wave vectors lie on a circle. 

In thin cells, where w/d, dqc ~ the three- 

dimensional amplitude equation (|l|) is not an accurate 
description anymore. But when solved numerically with 
B.C. (|) as a qualitative model, it predicts normal rolls 
in agreement with experiments. 

The theoretical analysis for thick cells seems to cor- 
respond well with the observations of Bohatsch and 
Stannarius in the "broad" intermediate frequency range 
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(within the high-frequency regime), where, near onset, 
the distribution of the wave- vector orientations is "broad, 
smeared out, and fluctuating" Unfortunately, no ex- 
perimental raw data for this range have been published, 
yet, which makes it hard to judge how well defined the 
maxima of these distributions arc, which they do never- 
theless find at small but finite obliqueness angles (« 10°), 
and which are difficult to understand from the theory 
presented here. 



III. BOUNDARY MODES WITH MODIFIED B.C. 

The experimental observation that, for frequencies 
above the intermediate "transition range" , the distribu- 
tion of the wave vectors near onset becomes narrower 
with maxima at obliqueness angles near ±45°, i.e. paral- 
lel to the surface anchoring, seems to be in contradiction 
with the theory based on the SM. 

A simple explanation, which should be examined ex- 
perimentally before coming to a final conclusion, would 
be the combination of two effects. The first effect is that, 
at higher frequencies, where the pattern wavelengths are 
shorter and wave-optical effects have to be taken into ac- 
count, the optical contrast of the shadowgraph method 
and the intensity of the reflexes in the optical far field be- 
come weaker. As a result, the apparent optical threshold 
of convection could lie slightly above the actual thresh- 
old. Wave-optical effects become important at pattern 
wavelength A of the order of magnitude of the geometri- 
cal average of the optical wavelength and the thickness 
of the convective layer = 0{w) = ©(Ad)-^/^ a condi- 
tion satisfied in the present case. The second effect is that 
the value of the reduced control parameter correspond- 
ing to the onset of nontrivial, nonlinear behavior involv- 
ing homogeneous soft modes scales like (A/d)^ or even 
(A/d)'^ |l7j, and thus decreases with increasing frequency. 
In combination of these two trends, the observed large 
obliqueness angles near threshold in thick cells at high 
frequencies might, in fact, correspond to already fully 
developed nonlinear effects, such as the chevron pattern. 
In particular the obliqueness angles above 45° (Ref [ pO[ , 
Fig. 11, 3500 Hz) would then be easier to understand. 

But one could also try to explain the observation of 
preferred obliqueness angles near ±45° as a result of some 
mechanism of wave- vector selection that removes the de- 
generacy of the problem. Mostly linear selection mech- 
anisms via perturbations from the boundaries shall be 
considered here. Nonlinear mechanisms will shortly be 
discussed in Section [V. 

If there is a linear selection mechanism for large 
obliqueness angles at threshold, at least one of the tree 
conditions (|^,^,^) must be invalid. The reason may be 
found in the truncation errors of the asymptotic expan- 
sion which leads to Eqs. the excluded effect of 
noise, or an insufficient description on the hydrodynamic 
level. All possible truncation errors become, by one way 



or another, large only for zq = 0{d). This will be demon- 
strated for one case below. If noise is the main cause for 
selecting zq (which is presumably the case at intermedi- 
ate frequencies), it is reasonable to assume zq to be ran- 
domly distributed over [0, d], leading again to zq = 0{d) 
on the average. Thus it can be concluded that, in fact, 
the description on the hydrodynamic level, in particular 
near the boundaries, is insufficient. 

For the eigenvalue problem considered above, the 
B.C. ^ turns out to be a somewhat singular case of 
the more general condition 



A — fidzA at z = —Zq 



(10) 



with a real parameter ^. To discuss the implications of 
/i ^ as one possible truncation error or as a possible 
outcome of an extended hydrodynamic description, the 
eigenvalue problem (||) with B.C. ([9HlO|) shall now be in- 
vestigated. This will also yield an analytic argument, 
robust with respect to perturbations, that supports the 
numerical observation that for ^ = the function e(zo) 
is monotonously decaying with zq. 

Up to normalization, there is, for fixed e, a single non- 
trivial solution A{z) = f{e, z) of Eq. (^ satisfying the 
upper B.C. (^). The lower B.C. ( |To| ) determines a dis- 
crete set of zo for given e or, respectively, e = e(zo) is 
determined implicitly by the lower B.C. as a multi- valued 
function of zqj i-e., by 



/(e(zo), -zo) = fif'{e{zo), ~zq), 



(11) 



where the prime on / denotes a differentiation with re- 
spect to the second argument. 

At a local minimum of s{zo) one has 9zj,£(zo) = 0, and 
differentiation of Eq. (pT|) with respect to zq yields 



/'(e(zo), -zo) = /i/"(e(zo), -zq). 



(12) 



By combining Eqs. ( pTp^ ) and the defining equation of 
/, Eq. (0), one arrives at 



2t2„2 .,2 



4d2 



(13) 



For nonzero /i the common factor / can be canceled since 
/ 7^ by Eq. (|l^). A simple algebraic relation between e 
and Zq (which determines the wave vector) for the critical 
mode is thus obtained. 

As /i ^ 0, two cases can be distin guis hed. In the first 
case Zq remains bounded. Then Eq. (y_3|) reduces asymp- 
totically to 



(14) 



Corresponding solutions exists only for ^ < 0. They can 
be approximated by /(e, z) ~ exp(z//i) and are strongly 
localized near the lower boundary. But due to the large 
values of |e| (as /.t — > 0) these solutions can not be ex- 
pected to be always physically meaningful. In the second 
case Zo increases over all bounds. Then small deviations 
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from eigenvalue and solution (J^) are sufficient to sat- 
isfy the B.C.. Thus e remains bounded and Eq. 
reduces asymptotically to 

In particular zq ^ oo as /i — > 0, supporting the numerical 
observation that e^zq) is strictly monotonously decaying 
for ^ = 0. 

It is easily verified that the analytic argument for 
— * oo as /i — > is, as far as physically interesting, ro- 
bust against additions of higher order corrections of the 
form g{(p'^)A to the r.h.s. of Eq. (||) or the replacement 
of the term by some more precise form h{(p^)d^A 

[assuming that h{ip'^) has no zeros for real ip, i.e., the 
PDE be nonsingular] . All other higher order corrections 
to Eq. (H) involve higher order derivatives (and additional 
boundary conditions) . If these corrections are small they 
are, in a first approximation, active only in a thin bound- 
ary layer and lead, effectively, to modified boundary con- 
ditions of the form (p^). Thus, taking the validity of the 
upper B.C. (||) for granted, a modified lower B.C. is al- 
ways involved in a localization near the lower boundary. 

To better understand the situation for fi =/= 0, Si numer- 
ical solution of the problem is required. Figure |l| shows 
numerical results for the values of s and zq at threshold 
[i.e. the minimum of £(20)] for finite fi. Although in the 
limit — > numerics break down, the asymptomatic re- 
sults are reproduced for small fi. The calculations also 
show that for fi > the minima of £(2:0) lie below the 
bulk (zq 00) value Ec as expected, but for < are, 
on the starred branch (s. Fig. |l]), above £c. A look at 
the function e{zo) (Fig. g) shows that for negative jj, the 
minima on the starred branch belong to the branch of 
e{zo) which approaches 3£c as zq — > cxo. It becomes also 
clear that the existence of the minima of £(2:0) for /i < 
is closely related to the existence of the thin boundary 
mode corresponding to the eigenvalue (|lj): The min- 
ima result from the avoided crossings of the spectrum 
of the modes localized near z = 0, which have decreas- 
ing eigenvalues for increasing zq, and the eigenvalue of 
the thin boundary mode, which increases with zq- If 
the thin boundary mode is an unphysical artifact of the 
boundary conditions the minima in the function e{zQ) 
are also unphysical. This point can be demonstrated by 
using the original B.C. (^) as the lower boundary condi- 
tion and adding, to model the mixed B.C. (|l0|), a term 
b{z + zo)dzA to the l.h.s. of Eq. (||), where b{z) is, e.g., 
given by 

^/ \ _ j 13.2£cW/i if z < 2.83w , 
^ ' 1 otherwise. ^ ' 

With these modifications the thin boundary mode — and 
therefore the avoided crossings and the minima of £(20) 
— are suppressed for /i < 0, but otherwise the spectra 
change only little (s. Fig. ||). 



Can higher order corrections to the boundary condi- 
tions that arise from the asymptotic expansion in the 
framework of the SM explain the preference of large 
obliqueness angles? From the figure in Ref. jl^ display- 
ing the hydrodynamic boundary layer of dielectric elec- 
troconvection, a value /i « —OAq~^ can be obtained. 
The sign follows directly from the fact that convection 
is suppressed by the no-slip boundary conditions for the 
hydrodynamic velocity field, the order of magnitude fol- 
lows from the geometry of the flow. By the sign of ^ a 
preference of large obliqueness angles (boundary insta- 
bility) is excluded. Even if the sign would change but 
the magnitude remain the same, Eq. ( p^ would give an 
unphysical value zq > d. 



IV. NONLINEAR SELECTION 

In principle a nonlinear selection of large obliqueness 
angles is also thinkable. At threshold, two kinds of non- 
linear mechanisms are active The first is the di- 
rect, local, nonlinear saturation of the pattern amplitude, 
the strength of which is measured by the coefhcient g in 
Ref. . This mechanism inhibits the coexistence of lin- 
ear modes of the form (^ if the overlap of the two modes 
along z is too large. Correspondingly, the wave vectors of 
coexisting linear modes can be expected to enclose some 
finite minimum angle. This could, for not too small w/d, 
lead to the preference of certain wave-vector combina- 
tions which fill the allowed range of « ±45° particularly 
well, but it does not imply that large obliqueness angles 
are generally preferred. 

The second nonlinear mechanism acts indirectly via 
the electric potential. At positions in the sample where 
\A\'^ is large, additional charge transport processes are ac- 
tive, which effectively increase (S'b in Ref. is positive) 
the electrical conductivity. Consequently, the strength of 
the electric driving field is reduced where \A\'^ is large and 
increased where \A\'^ is small, under the constraint that 
the total voltage drop is constant. But this effect is typi- 
cally weak compared to direct nonlinear suppression, and 
presumably does not change the qualitative properties of 
the planform selection problem. 



V. CONCLUSION 

It was shown that a twist of the nematic director - 
or another anisotropy axis in a different system - leads 
to a localization of the pattern in planes parallel to the 
director. When the twist is not too strong, the localized 
linear modes can be calculated in the framework of the 
amplitude formalism. In electroconvection, the z axis in 
Eq- (0) plays, apart form being the twist axis, another 
special role - the external electric field is parallel to z. 
But little would change if this additional anisotropy was 
missing (e.g. for reaction-diffusion patterns in uniaxially 
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anisotropic media). Only S^y and would then be equal 
by symmetry. 

For dielectric EC in twisted nematic (TN) cells the 
localization width w = 0{d/qc)^/'^ is given by Eq. (|^). 
When the material parameters of the liquid crystal are 
known, w can be calculated using the analytic approxi- 
mations for the coherence lengths and the wave number 
from Ref. |17|. As a result of the localization, the crit- 
ical wave vector is degenerate on a circular arc. The 
angle enclosed by the arc might be increased beyond 90° 
by using supertwist nematic cells (STN). One does then 
obtain a pattern-forming system with a geometry very 
similar to a two-dimensional isotropic system, but with 
a different topology (no invariance of the pattern under 
180° rotation) and very particular nonlinear interactions. 

The observation of normal rolls in thin cells and at low 
fre quen cies (large wavelength) by Bohatsch and Stannar- 
ius |2Cf| is compatible with — and presumably predictable 
from — the standard model. The transition to a state 
with a "broad, smeared out, and fluctuating" distribution 
of wave-vector orientations at higher frequencies can be 
understood from the wave-vector degeneracy. The the- 
ory also reproduces their observation that the transition 
takes place at some more or less fixed ratio \/d. This cor- 
responds to some fixed ratio w/d dX which the B.C. (^) 
become appropriate and the degeneracy sets in. Assum- 
ing that the observations at yet higher frequencies, where 
the wave vectors were better defined and aligned with the 
director anchoring, did not suffer from the technical dif- 
ficulties described in the beginning of Sec. [II, they seem 
to exhibit effects not described by the standard model. 
These effects seem to lead effectively to modified bound- 
ary conditions of the form ( p^ with /j, > 0. Neverthe- 
less this does not imply that, just as in the twisted cell, 
there are boundary modes also in the untwisted cell. So- 
lution of the linear problem for the untwisted cell with 
B.C. (|lO|), > 0, and a corresponding upper B.C. yields 
a critical mode covering the whole cell. Since the the- 
ory can explain the transition at the lower end of the 
intermediate-frequency range with fluctuating wave vec- 
tor but does not imply that there is another transition 
at an upper bound, it may be useful to investigate the 
two transitions independently also in experiments. The 
intermediate-frequency range is probably an independent 
regime, and not just as the crossover between the low- 
frequency and the high-frequency range. 

Conclusive information on the 3D structure of the 
high-frequency mode in untwisted cells might be obtained 
by investigating spectrally resolved shadow graphs or 
light scattering. 
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FIG. 1. (a) The threshold value of dielectric electroconvec- 
tion (relative to the bulk threshold in an untwisted nematic) 
in units of Sc [Eq. (^] and (b) zo (~ distance of the maxi- 
mum of the critical mode from the boundary) in units of the 
width w of the critical mode [Eq. (^)]. Both as functions of 
the parameter fj, in the boundary condition (^^. Stars and 
crosses are numerical solutions. The crosses correspond to 
a thin boundary mode. The dashed lines are the analytic 
approximations for small Eq. ( |l5| ) for (b). 




FIG. 2. The lowest eigenvalues e{zo) for (a) fj, = — 0.5w, 
(b) ^ = 0, and (c) ^ — 2.0 u;. The solid lines correspond 
to the eigenvalue problem Eqs. (^J^j^, the dotted lines to 
a modified model that does not contain the boundary mode 
leading to the large negative eigenvalues e{zo) for negative fi. 
For /I = solid and dotted lines are identical. 
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